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Abstract We study the steady-state structure of an 
accretion disc with a corona surrounding a central, ro- 
tating, magnetized star. We assume that the magneto- 
rotational instability is the dominant mechanism of an- 
gular momentum transport inside the disc and is re- 
sponsible for producing magnetic tubes above the disc. 
In our model, a fraction of the dissipated energy inside 
the disc is transported to the corona via these magnetic 
tubes. This energy exchange from the disc to the corona 
which depends on the disc physical properties is mod- 
ified because of the magnetic interaction between the 
stellar magnetic field and the accretion disc. According 
to our fully analytical solutions for such a system, the 
existence of a corona not only increases the surface den- 
sity but reduces the temperature of the accretion disc. 
Also, the presence of a corona enhances the ratio of gas 
pressure to the total pressure. Our solutions show that 
when the strength of the magnetic field of the central 
neutron star is large or the star is rotating fast enough, 
profiles of the physical variables of the disc significantly 
modify due to the existence of a corona. 
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1 Introduction 



stars: magnetic field 



Accretion discs are believed to be present in a wide 
variety of astronomical systems. The standard theory 
of geometrically thin accretion discs has provided suc- 
cessful framework for understanding the basic physical 
properties of these objects (Shakura & Sunyaev 1973). 
A solution for the structure of an accretion disc in an 
external magnetic field is desirable because accretion 
discs are frequently found around magnetic stars (Ver- 
bunt 1993; Warner 1995; Romanova et al 2003). The 
stellar magnetic fields interacts with the disc and can al- 
low the transfer of angular momentum between the disc 
and the star and this significantly affects the structure 
of the disc and the spin evolution of the star (e.g. Ghosh 
& Lamp 1978; Aly 1980; Wang 1987; Armitage, Clarke 
& Tout 1999; Schenker et al 2002; Dai & Li 2006). 
Matthews et al (2005) reformulated the standard disc 
solution to incorporate the effect of a torque from a stel- 
lar magnetic field. They studied the effects of varying 
the mass transfer rate, spin period and strength of the 
stellar magnetic field using their analytical solutions. 
Also, there are several simulations of the interactions 
between the stellar magnetic filed and accretion disc 
(e.g. Miller & Stone 1997; Romanova et al 2002; Ro- 
manova et al 2003). Miller & Stone (1997) investigated 
disc-star interaction for different geometries and stel- 
lar magnetic fields. Romanova et al (2002) described 
in detail the disc-star interaction for the case of slowly 
rotating stars and also cases of fast rotating stars. 

On the other hand, there are varying models of a 
hot corona above a cool disc (e.g. Haardtfe Maraschi 
1993). Begelman, McKee & Shields (1983) presented 
detailed calculations of the formation of the corona by 
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irradiation from the central source. Rozanska & Cz- 
erny (1996) and Rozanska et al (1999) investigated the 
effect of illumination by the central source on the ac- 
cretion disc for either a black hole as in AGN or a neu- 
tron star as in X-ray binary sources. It was found that 
the disc and the corona are separated by a boundary 
layer or heated upper region of the disc, i.e. the re- 
gion above the disc becomes stratified with increasing 
temperature and decreasing density (see also Rozan- 
ska et al 2002). Jimenez-Garate, Raymond & Liedahi 
(2002) calculated the vertical structure of an Accre- 
tion Disc Corona (ADC) and allowed the radius of the 
corona to extend to the same radius of the accretion 
disc. Low Mass X-ray Binaries (LMXBs) have two 
continuum components of radiation: simple blackbody 
emission from the neutron star plus comptonized emis- 
sion from an extended corona above accretion disc (e.g. 
Church & Balucinska-Church 1995). Church & Balu- 
cinska (2004) have demonstrated the importance of the 
size of the ADC to the correct description of Comp- 
tonization and they derived the comptonized spectrum 
of an LMXB based on the thermal Comptonization of 
seed photons. They estimated the radial extension of 
ADC in some LMXBs up to 10 9 cm. 

Generally, models for generating an ADC can be di- 
vided into two groups: (1) The dissipation of acoustic 
or magnetic energy flux generated by convective or vis- 
cous turbulence within the disc (e.g. Galeev, Rosner 
& Vaiana 1979; Merloni & Fabian 2002; Merloni 2003); 
or (2) the evaporation of disc material by the central 
source (e.g. Begelman, McKee & Shields 1983). 

Differential rotation of a standard accretion disc am- 
plifies the magnetic field. This may result in a magnetic 
active corona surrounding the disc (Burm 1986). These 
magnetic loops which thread the disc are analogous to 
the solar case. The problem of magnetic flux escape 
from an accretion disc has been studied by many au- 
thors (Stella & Rosner 1984; Hcyvaerts & Priest 1989; 
Shibata, Tajima & Matsumoto 1990; Romanova et al 
1998). 

The are reasons to believe that some parts of the 
ADC may find itself embedded in the field of the cen- 
tral object. Stollman & Kuperus (1988) assumed that 
the neutron stars in some of the LMXB have magnetic 
fields, which are capable of holding the accretion disc at 
some distance from the stars surface. They developed 
a model for the interaction of the magnetic field of the 
neutron star and the magnetic loops in the corona above 
the disc. Flaring interactions between accretion disc 
and neutron star magnetosphere has also been studied 
by many authors (e.g. Aly & Kuijpers 1990). Magnetic 
loops of the accretion disc which extend toward corona 
may contribute to angular momentum transport in the 



disc (Burm & Kuperus 1988). However, the possible 
effects of the corona on the structure of a thin accre- 
tion disc under the influence of a stellar magnetic field 
have not been studied to our knowledge. In this paper, 
we extended Matthews et al (2005) analysis by con- 
sidering the corona of the disc as well. The following 
question is addressed: what is the effect of the corona 
on the steady state structure of a magnetically torqued 
thin accretion disc? We obtain a set of analytical solu- 
tions for a thin accretion disc with a corona, to which 
a magnetic torque due to the central object (e.g. neu- 
tron star) is applied. We perform a parameter study to 
explore how physical variables of the disc such as den- 
sity or temperature may be modified due to changes 
in stellar magnetic field or stellar spin. In section 2, 
we discuss the basic assumptions and present the main 
equations of the model. Typical physical properties of 
the solutions will be studied in section 3. In the final 
section we conclude the paper by a summary of our 
results. 



2 General formulation 

We assume that all angular momentum transport takes 
place in the disc and the mass accretion rate M is 
constant with radius and time. However, the micro- 
physics mechanism of the angular momentum transfer 
remains unknown. Shakura & Sunyaev (1973) replaced 
all the missing physics by a parameter a. This ap- 
proach has been widely used for studying the dynamics 
and structure of accretion flow. A promising mecha- 
nism for driving the turbulence responsible for angular 
momentum and energy transport is the action of the 
magneto-rotational instability (MRI) that is expected 
to take place in such discs (Balbus & Hawley 1991). 
But for a disc-corona system, it is typically thought 
that the viscous stress, assumed to be magnetic in na- 
ture, transports angular momentum and initially ran- 
domise the gravitational binding energy near the mid- 
plane. The magnetized fluid elements, which are buoy- 
ant with respect to their surroundings, dissipate above 
the disc. Merloni & Fabian (2002) and Merloni (2003) 
presented very detailed discussions about dissipation 
in the corona and its relation with angular momentum 
transport in the disc itself. Our disc-corona model is 
developed along the line proposed by Merloni & Fabian 
(2002) and Merloni (2003). 

We consider a more general prescription for the vis- 
cous stresses r x( j, (Taam & Lin 1984; Watarai & Mi- 
neshige 2003; Merloni & Nayakshin 2006): 

= -aop'-^p^, (1) 
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where ao and < \i < 2 are constants. Also, p is the 
sum of the gas and radiation pressures. Phcnomcnolog- 
ical models generally assume that at each radius, a frac- 
tion / of accretion energy is released in the reconnected 
magnetic corona. Assuming that turbulence inside the 
disc is MRI-driven, such a fraction / of the binding 
energy is transported from large to small depths by 
Poynting flux, Merloni & Nayakshin (2006) estimated 
the fraction / as 

/ = v/2«o^/ 2 , (2) 

where (3 is the ratio of gas pressure to the total pres- 
sure. So, in this model, this fraction / is not a free 
parameter. Actually, it is very difficult to present a de- 
tailed analytical model which gives us the fraction / as 
a function of the disc physical variables structure and 
MR1 (even without central magnetic torque). We ap- 
preciate on possible effect of central magnetic torque on 
energy exchange /. But it is beyond our study because 
of severe limitations of our phenomenological approach. 
But we would like to emphasize that the amount of the 
energy exchange corresponding to solutions with cen- 
tral magnetic torque (our solutions) and without cen- 
tral torque is not the same. Because for solutions with 
a central magnetic torque, we have an extra magnetic 
torque in angular momentum equation (see equation 
(J7J) which dramatically changes the solutions compar- 
ing to the solutions without that term (see solutions 
of Matthews et al. (2005)). In fact, equation ([2]) pre- 
scription is based on MRI inside of the disc, irrespective 
of the (non)existence of external magnetic torque. So, 
with central magnetic torque its effect on the energy ex- 
change appears by changing disc structure (through an 
extra magnetic term in angular momentum equation). 

The rotation curve is dominated by a Newtonian 
point mass M, as relativistic effects are only important 
at small radii. Thus, the rotational angular velocity of 
the disc is Keplerian, i.e. f^K = ^GM/R 3 . Having 
equation (fTJ) as a prescription for the viscous stresses 
and equation ^fy as an expression for the fraction of 
power dissipated in the corona, we can write the basic 
equations of the disc. Vertical hydrostatic equilibrium 
of the disc implies 

P _ & K H 

S~ 2 • [6) 

In order to write the azimuthal component of equa- 
tion of motion, the specific magnetic torque of the cen- 
tral neutron star can be written as 

A=^ V2 , (4) 



where tA is the time-scale on which the local disc mate- 
rial gains angular momentum. This torque appears be- 
cause of the interaction of the field lines of the central, 
rotating magnetized neutron star with its surrounding 
partially ionized inner accretion disc. Actually nei- 
ther the physics of angular momentum transport by 
the magnetic field of the central star, nor the mecha- 
nism of penetration of the magnetic field into the ac- 
cretion disc are well understood. Nevertheless we can 
present some physical prescriptions for the magnetic 
torque of the central neutron star. Actually, there are 
some interactions between the magnetic tubes above 
the disc and the magnetic field of the neutron star. 
However, we concentrate on the steady structure of the 
disc and such complicated interactions are neglected in 
our model. Formation of the corona is a direct conse- 
quence of the magnetic tubes of the accretion disc. On 
the other hand, the physical properties of the accretion 
disc are modified because of the existence of the mag- 
netic torque and the dissipated energy into the corona 
depends on the structure of the disc. Thus, the effect of 
the stellar magnetic field on the magnetic tubes is in- 
cluded via its torque on the disc indirectly in our model. 
Figure [T] show schematically the magnetic field of the 
neutron star and the magnetic tubes of the disc. The 
magnetic structures in the corona of the disc should be 
of a small scale, because turbulent eddies in the disc 
are believed to have the disc thickness as a character- 
istic size, as should the magnetic loops which pop out 
as a result of magnetic buoyancy. However, there is a 
possibility of forming larger magnetic loops above the 
disc by reconnecting small scale loops (e.g., Heyvaerts 
& Priest 1984). Analysis of such complexities is beyond 
this paper. 

In equation (jl]), the torque time-scale tA is equiva- 
lent to the magnetic time-scale i ma g- Magnetic field of 




Fig. 1 Schematic of the magnetic field of the neutron star 
and magnetic tubes of the accretion disc. The truncation 
and the corotation radii are marked by Rt and R co , respec- 
tively. 
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the neutron star has a dipolar form, i.e. B z ~ |/i|i? -3 
where [i is the magnetic moment of the star. Depend- 
ing on the magnetic interaction model, the magnetic 
time-scale can be parametrized in the form (Livio & 
Pringlc 1992; Campbell & Hcptinstall 1998; Matthews 
et al 2005) 



2E 



#7+2 



Pm [{R/Rcof 2 - 1] ' 



(5) 



where 7 is determined by the magnetic interaction 
model (e.g., Campbell & Heptinstall 1998) and (5m is 
calculated by the mass and the magnetic moment of the 
star, 



M 



l/'-l 



2nVGM 
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In our fully magnetized case we have 7 = 7/2 (Livio 
& Pringle 1992; Matthews et al 2005). Also, in the 
above equation, the corotation radius is represented by 
R co which is by definition the radius of the accretion 
disc at which, the star and the disc rotate at the same 
angular rate (Figure [T]). Note that the dipole and spin 
axes are taken to be parallel to z, and the whole system 
is assumed to be axisymmetric. 

Now, we can write the azimuthal component of the 
equation of motion as 



dR 



RdR 



(7) 



where wr and v v are the radial and the rotational ve- 
locities of the disc. We note that the disc is rotating 
with Keplerian profile, i.e. v v = vk- The r — (p compo- 
nent of viscous stress tensor is prescribed by equation 
(TTJ). Also, the parameter A is representing the specific 
torque of the central neutron star which is given by 
equation After mathematical manipulations, we 

can write equation jJJ as 



7-0/2 



3Cl K Mq 4 v, 
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where, 
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M i?T( 7 -2) g 4 ' 
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The truncation radius Rt denotes the inner boundary 
of the accretion disc. A common feature of all magnetic 
star-disc interaction models is that the stellar magneto- 
sphere is strong enough to disrupt the accretion disc at 



some radius Rt above the surface of the star. There are 
several methods to locate the truncation radius. We use 
the approach of Matthews et al (2005). They presented 
an approximate relation for 7 = 7/2 as 



R 



Q 



(9) 



where the nondimcnsional parameter Q is determined 
by the ratio /3m, accretion rate and the corotation ra- 
dius, i.e. 



2nR^(^f). 
M 



(10) 



The truncation and corotation radii are shown schemat- 
ically in Figure [TJ 

The energy equations is given by 



(11) 



where as we discussed / = ^/2ao/3 M / 2 . Also, we can as- 
sume that the vertical transport of heat is by radiative 
diffusion which implies the midplane and the surface 
temperatures are related by 



r = (^s) 1 /4 Tcff , 



(12) 



where k is the opacity coefficient. In order to limit the 
parameter space in our analysis, we simply assume the 

electron opacity for this coefficient, i.e. k = K e = 0.4 

2 —i 
cm g . 

Now, equations ©, {5J and (jTTJ) are the main equa- 
tions which enable us to find p and T as functions of R 
and (3 and the other input parameters. Thus, 



167r 2 agc/cB 



4crQ K 1/2 



>-l/3 



(1 - 0)V3 P 



4<rn K ^_ 1/2( 167T 2 o; 2 cfcB 



(13) 



3a/j, m mB_M 2 fl^q s v 2 



(1- A /W^)V 2 fa _ BwVia 
(l-/3) 2 /3 P 

There is an algebraic equation for (3 as follows 



(14) 



fc B ,4crfi KN _ 3/2/ 16TT 2 alck B - 5/3 



^ m m H 3/sao 



-r 3/2 ( 
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Fig. 2 Profiles of the physical variables for a disc with a corona (solid lines) and without a corona (dashed lines) vs. radial 
location in the disc with fi = 1.0, /i m = 0.6, k = 1, Mi = 1.4, Mis = 1, P Bp in = 1 s and B s = 10 12 G. Each curve is labeled 
by its corresponding viscosity coefficient ao- Also, in the online version of the paper, curves are shown by different colors 
depending on the viscosity coefficient. The top left-hand and right-hand plots show the surface density and the temperature 
of the accretion disc, respectively. The ratio of the thickness of the disc and the radius is shown in the bottom left-hand 
plot. Finally, the bottom right-hand plot shows profile of the ratio f3 of gas pressure to total pressure. Generally, we see 
that the corona becomes more effective when the viscosity coefficient ao has larger values. 




Fig. 3 Profiles of the physical variables for a disc without a corona (dashed lines) and with a corona (solid lines) vs. radial 
location in the disc with cto — 0.3, /u m = 0.6, k = 1, Mi = 1.4, Mis = 1, P sp in = 1 s and B s — 10 12 G. These plots show 
dependence of the physical variables to the viscosity exponent /j,. Typical behaviors of the ratio H/R and the ratio of the 
gas pressure to the total pressure, /3, for the case without a corona are not very sensitive to the exact value of the viscosity 
exponent. But in the case of nonexistence of a corona, profiles of the surface density and the temperature of the accretion 
disc are not independent of the parameter [i. As this exponent increases, both the surface density and the temperature 
increase. 
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Fig. 4 Profiles of the physical variables for a disc without a corona (dashed lines) and with a corona (solid lines) vs. radial 
location in the disc with ao = 0.3, /i = 1.0, /U m = 0.6, k = 1, Mi = 1.4, Mi6 = 1 and P sp m = 1 s. Effects of the magnetic 
field of the central neutron star are shown in these plots. Each curve is labeled by the magnetic field at the surface of the 
neutron star, i.e. B s . As the magnetic field increases, the possible differences between the solutions with a corona and 
without a corona become more significant. The bottom right-hand plot shows profile of the ratio /3 of gas pressure to total 
pressure. 




Fig. 5 Profiles of the physical variables for a disc without a corona (dashed lines) and with a corona (solid lines) vs. radial 
location in the disc with Qo = 0.3, fi = 1.0, ^ m = 0.6, k = 1, Mi = 1.4, Alia = 1 and B s = 10 12 G. The bottom right-hand 
plot shows profile of the ratio j3 of gas pressure to total pressure. The central magnetic neutron star is rotating at different 
speed and each curve is labeled by the rotational period of the star. We see that the truncation radius increases as the 
central neutron star is rotating slower. 
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where /x m is the mean particle mass in units of the hy- 
drogen atom mass, mjj- The other constants have their 
usual meanings. 

In order to study the behavior of our solutions, it is 
more convenient to introduce dimensionless variables. 
For the central mass M, we introduce Mi = M/Mq 
and for opacity we assume k = k/k £ . The radial dis- 
tance from the star, in the plane of the disc, is rep- 
resented by R in units of 10 10 cm. So, wc introduce 
no = i?/(10 10 cm). Also Mie = M/(10 16 gs _1 ) repre- 
sents the mass accretion rate in units of 10 16 g s 

We can rewrite our solutions as 



p = 2.63 x 10- 10 ao 1/2 k 3 ^M^Mf 6 q 8 v 2 r^ 15/i 



x/3 (8-M)/4 (1 _ ^-2^ _ v /2a /3^/2)3/2 5 
p = 40.47 X a - 5/6 / t 1 /2 Ml 1 3/ 12 M l 4 /3 g 16/3 w 4/ 3ri -13/4 



(16) 



x/3 (8-5 M )/i2 (1 _ /3) -2/3 (1 _ v / 2 a ^ /2 ) 1/2 , (17) 
I = 0.24 x a^K-^M^M^q-^v- 1 ' 3 



x/3 -(8 +M )/12 (1 _ p) 2/3 {1 _ v / 2ao ^/2)-V2 j (18) 

and the ratio (3 is obtained from the non-dimensional 
form of equation (fT5)) , that is, 

4.69 x 10- 7 al /6 k 3 / 2 Ml /12 Mtl 3 q w/3 v A / 3 r^ /4 



x^ 8+fi)/12 (l - /3) 5/3 (l - y / 2a p^ 2 ) 3/2 = 1. (19) 
We can also calculate the surface density as 



£ = 2.12 x 10- V /3 kMi 17/12 M^ 2 °/V/Vo 11/4 



x/3 (4-a0/3 (1 _ ^)-4/3 (1 _ ^/2a ^/ 2 ) 1/2 



(20) 



Typical disc structure will now be illustrated, using 
the full disc solutions quoted in equations (flT>]) . (TIT)) and 
(|18|) . However, the physical variables depend not only 
on the radial distance but also on the ratio of the gas 
pressure to the total pressure which can be calculated 
at each radius from algebraic equation (|19|) . Having 
analytical solutions, we can study physical properties 
of our disc and corona system. In the next section, wc 
will do a parameter study of our solutions. 



3 Analysis 

We fix the central mass, the mass transfer rate, the 
opacity and the mean molecular weight, respectively 
as Mi = 1.4, Mi6 = 1, k = 1 and fx m = 0.6. The 
other input parameters are changed to illustrate their 
possible effects on the physical properties of the sys- 
tem. Typical behaviors of the solutions for disc with 
a corona (solid curves) and without a corona (dashed 
curves) are shown in all subsequent figures. Each curve 
is labeled by its corresponding parameter, though color 
versions are available in the online version of the pa- 
per. Note that the radial extension of the plots are up 
to 10 10 cm, significant effects of the corona and the cen- 
tral magnetic field are seen at radii less than 10 9 cm, 
generally. 

Figure [2] shows typical behavior of the physical quan- 
tities of the accretion disc when the viscosity coefficient 
ao changes from low value 0.03 to high value 0.3. The 
profile of the surface density is shown in the top left- 
hand plot of Figure^ When the viscosity coefficient in- 
creases, the surface density decreases both in the cases 
with and without a corona. This behavior can be ex- 
plained simply by the fact the the accretion rate is con- 
stant in our model. So, the radial velocity (which is 
proportion to viscosity coefficient) is inversely propor- 
tional to the surface density. Thus, when the viscosity 
coefficient increases, the radial velocity increases too. 
But the surface density decreases with the viscosity co- 
efficient. However, the corona significantly increases 
the surface density. The top right-hand plot of Fig- 
ure [5] shows the presence of a corona serves to cool the 
disc. In fact, corona extracts some of the dissipated 
energy inside the disc. When there is not a corona, 
all of the dissipated energy heats the disc. But with 
a corona some of the dissipated energy transfers from 
the disc to the corona and so, the disc becomes cooler. 
Generally, when the viscosity coefficient increases the 
disc temperature decreases both with and without a 
corona. In fact, the fraction of the dissipated energy 
into the corona is directly proportional to ao according 
to equation @. Thus, as the profile of the tempera- 
ture shows, the existence of a corona leads to a cooler 
disc and the corona becomes more effective when the 
viscosity has larger values. The bottom left-hand plot 
of Figure [2] shows the ratio of the disc thickness to the 
radius. The corona causes the disc to become thinner 
than the case without a corona. It justifies the thin disc 
approximation even in the presence of a corona. The 
ratio of the disc thickness to the radius decreases when 
the viscosity coefficient increases. Hydrostatic equilib- 
rium of the disc in the vertical direction implies that the 
thickness of the disc is in proportion to the sound speed 
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(i.e. temperature). Since the existence of the corona re- 
duces the temperature of the disc, larger values of the 
viscosity coefficient imply a thinner disc. 

The bottom right-hand plot of [2] shows the ratio (3 
of gas pressure to the total pressure of the disc. In the 
absence of a corona, when the viscosity coefficient in- 
creases, the ratio of gas pressure to the total pressure 
decreases (i.e., radiation pressure increases). However, 
in a disc with a corona, the parameter (3 increases as 
the viscosity coefficient increases (i.e., radiation pres- 
sure decreases). Generally, the presence of a corona in- 
creases the ratio of the gas pressure to the total pressure 
in comparison to the same disc but without a corona. 
Clearly, the ratio f3 tends to unity at large radial dis- 
tances. On the other hand, the fraction / of energy 
exchange is directly proportional to the ratio (3 (see 
equation Thus, the fraction / reaches to its max- 
imum value when the ratio (3 tends to the unity. Thus, 
the maximum value of / is \/2ao according to equation 

Figure [3J shows the typical behavior of the physical 
quantities of the disc for different values of the viscosity 
exponent /x. The top left-hand plot of Figure [3J shows 
that the surface density of the disc increases when the 
viscosity exponent increases in both cases with a corona 
and without a corona. The presence of a corona en- 
hances the surface density for all of the input parame- 
ters of the viscosity exponent. However, in the case of 
a system without a corona, the profiles arc more sensi- 
tive to the variations of the viscosity exponent /z. The 
behavior of the temperature is shown in the top right- 
hand plot of Figure [3J The temperature of the disc in- 
creases when viscosity exponent increases for both cases 
with and without a corona. Generally, the presence of 
a corona serves to cool the disc, but as the viscosity 
exponent increases the disc becomes warmer in partic- 
ular in the inner parts. The bottom left-hand side of 
Figure [3] shows that the ratio of disc thickness to the 
radius increases with increasing values of the parame- 
ter n, though the enhancement is not significant. The 
presence of a corona causes the disc to become thin- 
ner in comparison with the same disc without a corona 
for all of the input parameters of the viscosity expo- 
nent. The bottom right-hand plot of Figure [3J shows 
that the corona causes (3 to increase in comparison to 
the case without a corona. When the viscosity exponent 
increases then the ratio (3 increases for a disc without 
a corona but decreases for the case of a disc with a 
corona. 

The measured values of the magnetic strength of the 
neutron stars are in the range 10 8 5 G to 10 12 G (Dewey 
et al 1986; Frank, King & Raine 2002). This implies 
that the magnetic fields of the neutron stars in the 



LMXB are at least as strong. It was shown that the 
magnetic field of neutron stars probably do not decay 
below 10 9 G or if they do the timescale is larger than 
10 9 yr (Van den Heuvel et al 1986). Typical behav- 
iors of the physical variables with different strengths of 
the stellar magnetic field are shown in Figure [U The 
top left-hand plot of Figure H] shows the surface density 
of the disc. The surface density of the disc increases 
when the magnetic field increases for both cases with a 
corona and without a corona. The presence of a corona 
enhances the surface density. The top right-hand plot 
of Figure Q] shows the disc temperature increases when 
magnetic field increases for both cases with a corona 
and without a corona. Presence of the corona decreases 
temperature of the disc for all of the input parameters. 
Generally, we see that the corona becomes more effec- 
tive when the magnetic field has larger values. The 
bottom left-hand plot of Figure 0] shows that the ratio 
of disc thickness to the radius increases when magnetic 
field increases. The corona causes the disc to become 
thinner in comparison with the same disc without a 
corona, also when magnetic field is stronger, the corona 
becomes more effective. The bottom right-hand plot of 
Figure H] shows the ratio of gas pressure to the total 
pressure of the disc. The ratio (3 decreases when the 
magnetic field increases for both cases with a corona 
and without a corona. The corona causes ratio of gas 
pressure to the total pressure of the disc to increase for 
all of the input parameters. 

The effect of the spin period on the physical param- 
eters is displayed in Figure [5] The profile of the surface 
density is shown in the top left-hand plot of this figure. 
The surface density of the disc decreases when the spin 
period of the star increases for both cases with a corona 
and without a corona. For all of the input parameters, 
the presence of a corona enhances the surfaces density. 
The right-hand plot of Figure [5] shows the corona de- 
creases temperature of the disc in comparison with the 
same disc without a corona, irrespective of the spin pe- 
riod. Also, the temperature of the disc decreases when 
the spin period increases and clearly this reduction of 
temperature is more evident when the star is rotating 
slower. The bottom left-hand side of Figure [U shows 
that when spin period increases then the ratio of disc 
thickness to the radius decreases for all cases. In this 
regards, the corona becomes more effective when the 
star is rotating faster. Also the presence of the corona 
causes the disc becomes thinner for all of the input 
parameters. The bottom right-hand plot of Figure \5\ 
shows that the presence of a corona enhances the ratio 
of gas pressure to the total pressure of the disc. The 
ratio of gas pressure to the total pressure of the disc 
increases when spin period increases and even it tends 
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to unity for spin periods longer than 5 seconds which 
implies negligible radiation pressure. Generally, we see 
that the corona becomes more effective when the spin 
period has smaller values. 

4 Conclusion 

We obtained analytic solutions for the thin accretion 
disc with corona, to which a magnetic torque due to 
the central object (e.g. neutron star) is applied. We 
assumed that all angular momentum transport takes 
place in the disc and the mass accretion rate M is con- 
stant with radius and time. In our simple model there 
is only energy exchange between the corona and the 
underlying optically thick disc. We showed that the 
existence of corona affects structure of the disc which is 
under the influence of the central magnetic field. Cen- 
tral in our study is the assumed existence of magnetic 
loops anchored in the disc and extending into a corona 
on both sides of the disc. Such loops are believed to 
arise by MRI inside the disc. We showed that the pres- 
ence of a corona serves to cool the disc because of the 
energy transported from the disc to the corona. This 
reduction of the temperature causes the disc becomes 
thinner. The surface density increases when we con- 
sider the corona in comparison to the same disc with- 
out a corona. Our solutions indicate that the presence 
of corona increases ratio of the gas pressure to the total 
pressure. 

Although we considered a generalized prescription 
for the viscosity, profiles of the physical variables are 
weakly depending on the viscosity exponent in both 
cases with a corona and without a corona. But behav- 
iors of the solutions are highly depend on the viscosity 
coefficient. The effect of corona is more evident as vis- 
cosity tends to high values. Note that we can not di- 
rectly compare our solutions with Matthews et al (2005) 
study. In fact, our viscosity prescription and viscosity 
coefficient are not similar to the analysis of Matthews et 
al (2005). Moreover, we considered both radiation and 
gas pressures. Another important factor in our model is 
the strength of the central magnetic field at the surface 
of the neutron star. For low level strength of the sur- 
face magnetic field, the effect of corona on profiles of the 
system is negligible. Also, in this case, the disc is in gas 
pressure dominated regime. Corona significantly mod- 
ifies physical profiles of the system with increasing the 
strength of the central magnetic field. The spin period 
of the central neutron star is also another important 
factor in our model. As the star rotates faster, the ef- 
fect of corona on the behaviors of the physical variables 
becomes more significant. These findings clearly show 



the importance of corona in LMXBs which can not be 
neglected in modeling of such systems. We note that 
the results are obtained based on only energy exchange 
from the disc to the corona. Further improvements of 
our model will include a more accurate description of 
disc-corona interactions, such as mass exchange or even 
angular momentum transport. Also, it will be interest- 
ing to calculate modifications to the disc spectrum due 
to the existence of corona according to our analytical 
solutions. 
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